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Abstract 

Model computations show that changes of sampling interval introduce only OJ 1 cm ^changes 
whereas zero padding provides an improvement of more than 5 cm in the FFT generated g 
For the GPS survey of Franklin County, Ohio, the parameters selected as a resu t of mo 
plu„., .11- b, reaction in loci data ",«».».«•. -«• th « "» 

accuracy when tapering and padding is applied. The results are shown in tables. 


Introduction 

rpi f .,. win(7 is a br ief description of computational modeling carried out in order to obtain optimal 

results from the use of the Fast Fourier Transform (FFT) technique for local geoid computation. These 
. j id £ j a. moat favorable parameters for local geoid computation using gravity 

rr:x n ^ — ■ b °* k r rr ,nd •tsit 

in these studies provides the opportunity to test some interesting aspects of the FFT technique. 


Model Description 

Fo, detad. ... Nagy{1988). The 

analytical expressions for the potential, U, and the gravity, Ag, for a single prism are g.ven below 
(Nagy, 1980) : 


U 


kp 


x y ln(z + r) + y«ln(i + r) +zxln(y + r) 


--x 2 arctan V — - \y 2 arctan ^ - \z 2 arctan ^ 

2 xr2 V r L 1 J 


A g = kp 


where 


i ln(y + r) + y ln(i + r) - z arctan 


xy] 

zr\ 


• = x/x 2 + y 2 + z 2 


The negative ot the potent,,, divided by the 

Summing up the required quantities or a pru, compared The difference in clearly the 

which the result of the various numerical computations can be c P metho d 

error ot the numerical procedure. In thi. c«, the error generated by the FFT method. 


Geological Survey of Canada Contribution No. 34788 


92 



NAGY AND Fury : FFT — Local Gravimetric Geoid Computation 


Effect of Sampling Interval 

As the transfer function, i.e. the function used to weigh the gravity anomalies to produce the geoidal 
height, is relatively flat as compared, for example, with the functions used in calculating the deflections 
of the vertical, or the vertical derivatives, one expects no large changes associated with the changes in 
the sampling interval This has been confirmed with model computations. Different sampling intervals 
between 1 and 15 km covering the same area produced only 0.3 cm change in geoidal height. For this 
reason, the sampling interval does not seem to be of major concern in local geoid determinations. 


Effect of Padding 

The Fourier method assumes periodicity, i.e. the field given in a two-dimensional array is repeated in the 
frequency domain around the central part in both dimensions and introduces the so called leakage into 
the computations, causing unwanted errors. To partially compensate for this error, the technique known 
as padding is used. Padding consists of putting zeros around the values of the input matrix, practically 
doubling the dimensions. For the model using 5 km sampling interval and a 26x26 grid, the gravity was 
practically zero at all boundries. The model geoid over this grid has a span (difference between maximum 
and minimum) of 74.9 cm. The use of FFT on the corresponding gravity anomaly produced a geoid with 
the span of 67.3 cm i.e. an error of 7.6 cm. Carrying out the zero padding to generate a matrix of 50x50 
resulted in a different geoid with a span of 73 cm. This means that doing only zero padding, the error 
was reduced from 7.6 cm to 1.9 cm. This is a far greater change then produced by varying the sampling 
interval. Here the great importance of modeling is stressed. The results of computations without and 
with padding are different. However without the knowledge of the exact model values^ one would not be 
able to draw any conclusions. In the case of modeling, the comparison with the exact values makes it 
obvious which computation gives the better result. 


Effect of Tapering 


Normally the gravity anomalies at the boundaries are not zero. In order to have a smoother transition, the 
techique known tapering is used. The purpose of tapering is to bring down the non-zero gravity values at 
the boundary smoothly to zero. There are various ways of achieving this, but model computations show 
that the particular method used for tapering is not critical. Table 1 summarizes some numerical results 
with various combinations of tapering with zero padding. The input matrix was generated at a 1 km inter- 
val, consisting of an array of 62x62 (used as reference), covering the central part of the anomaly field, with 
reasonably large non-zero values at the boundaries. All geoidal height related quantities are given in cm. 

Table 1 


Array 

Size 

Tapering 

% 

Padding 

% 

Geoidal Heights 

Residual Errors ] 

min 

max 

Span 

RMS 

62 

0 

0 

-24.55 

23.89 

40.73 

7.43 

70 

6 

0 

-22.31 

26.74 

27.30 

4.87 

90 

6 

16 

-15.98 

35.17 

22.51 

2.96 

90 

22 

0 

21.22 

34.68 

11.98 

2.54 

no 

22 

16 

-17.38 

42.76 

8.54 

1.77 

110 

38 

0 

-21.71 

42.34 

16.63 

2.96 

130 

38 

17 

-17.13 

49.69 

10.97 

2.16 

110 

6 

32 

-12.22 

41.78 

20.97 

2.66 

130 

22 

33 

-13.17 

48.99 

7.76 

1.64 
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Numerical Errors 

„ is well known th.t the competed vulu.s toward the horde, of the ^me e^neour Modehng 

r.rrs.ndSr^ 

top in residual errors uft„ reducing the urr.y .»« by about 10%, thereafter no sign, Scant reduction 
errors occur. Obviously, the results are again model dependent. 

Practical Application 

Based upon the results of model computations, the FFT method has been used to calculate relative 
geoidal heights for the Franklin County GPS survey. The calculation was done in two steps . 

. the regional component was calculated from the OSU86F truncated to n = 36, 

. the local component was derived by applying the FFT technique to the residual gravity held, wh, <* 
was gridded at 5' intervals resulting in an array of 192x192 providing the desrred coverage 

area of interest. 

The guoidul height difference » the sum of the globul geoid nnd column [1]. This vine will be used Inter 

Computation, are shown in Table 2. Column |.| is the dime, application of 
FFT all other solution, are the change, with respect to this solution. The dramatic reduction of error, 
by L combLd affect of tapering and padding (for enample, solution |8] vs. solution |4|) „ readily 

recognizable from Table 2. 

Table 2 



Base line 


Rhodes — 


18-83 — 


18-83 
Britton - 


Hoover - 


18-83 -► 


Jackson * 


Smith — ► 


Smith — ► 


Legend : 


■+ Clark 


Clark 


Rhodes 
-» 18-83 


-► Clark 


Shannahan 


Britton 


Jackson 


Hoover 


Global 

geoid 


-8 


Residual FFT Geoids 



[21 

[31 

[4] 

[5] 

[61 

| 7 ] 18 ] 

-.3 

.5 

2.4 

1.1 

.4 

.6 -.4 

-.4 

-.1 

1.6 

.9 

.3 

.2 -5 

-.1 

-.6 

-.8 

-.2 

-.1 

-.4 -.1 

-.1 

-.8 

-2.9 n 

-1.0 

-.1 

-.7 .7 

.5 

— 

.1 

-.8 

-1.1 

-1.1 

-.3 -5 

-.6 

Li 

5.0 

2.7 

.5 

.8 -1.3 

-.3 

1.9 

7.2 

2.9 

1.1 

2.2 -.7 

.4 

.9 

2.4 

.5 

-.5 

.4 -1.2 

-.9 

1.8 

9.1 

4.4 

1.9 

2.4 -1.2 


Array Rem arks — 

192 8.0° border around baselines 

144 5.5° border around baselines 

72 3.0° border around baselines 

52 20 rows and columns removed at south and 

east (simulating lack of data at shore lines) 
72 zero padding for removed data 
144 50% zero padding on [5] all around 

72 20% tapering on [4] all around 

144 50% zero padding on [7] 
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The relative geoidal height computations were repeated next by truncating the OSU86F global model 
to n — 180 and both results (i.e. the n — 36 and n = 180) were then compared to values derived from 
GPS surveys and leveling (Table 3). The GPS survey used in these comparisons was reported earlier 
by Engelis et al. (1984) and Kearsley (1985), and are listed as OSU and AUS respectively in Table 3. 
Relative geoid heights (GPS) were derived at the National Geodetic Survey of USA (Fury,1985); the old 
values were given in the above cited references. 


Table 3 


Base line 

Length 

Relative Geoid Heights j 

GPS 

old 

OSU 

AUS 

36 

180 

Rhodes — ► Clark 

10 

-5 

(-7) 

-3 

-6 

-5 

-3 

18-83 Clark 

11 

-18 

(-19) 

-14 

-19 

-16 

-14 

18-83 — ► Rhodes 

4 

-13 

(-13) 

-11 

-13 

-12 

-12 

Britton -♦ 18-83 

13 

-21 

(-19) 

-21 

-19 

-21 

-22 

Hoover — ► Clark 

10 

19 

(19) 

12 

15 

12 

12 

18-83 — ► Shannahan 

22 

-18 

(-25) 

-11 

-19 

-13 

-10 

Jackson — ♦ Britton 

24 

-4 

(1) 

0 

-5 

-2 

0 

Smith — + Jackson 

14 

58 

(32) 

50 

63 

51 

52 

Smith — * Hoover 

35 

4 

rfis) 

3 

5 

0 

4 


Conclusions 

The numerical experiments presented here confirm the effectiveness of the FFT method for local gravi- 
metric computation and show some of the results which can be obtained from model computations for 
use as guidelines in practical applications to obtain the best result from the FFT technique. 
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Abstract 

Gravimetric geoids have been computed for the °*^ a 18 o and 

terranean Sea between latitude 32 and 3b a h y ross _ over ad _ 

22° using EFT c °^“ n : n ? \£e area showed for both gravi- 

justed SEASAT and GEOSAT ~ * differences to be 

Metric geoids the deviation 

0.20 m and 0.15 m, respec Y the collocation geoid 

of the difference between ^respectively ■ This quite large 

heights were 0-82 a • differen £ data sampling and noise 

difference may be due to th fully explain- 

weighting used by the two methods, but is not yet Y 

ed. 

1. introduction 

.u i„ i QQO ' ties the ERS-1 and the Topex/Poseidon sa- 

tellies will be launched, both equipped with a radar altimeter. 

^usefulness of bhe altimeter data or oceanogra^hrc purpo^ 

will be greatly improved, if ve are a ' ^ investigated, i.e. a 

height reference surface for ^ rnean 9 that height differen- 

regiona 1 , relative, geoi . y n the va lues may be affected 

ces are precisely ‘b better, if we could com- 

pute ^Tbso^rte-gioid^b^ this will require that e.g. a global 
gradiometric satellite mission is carried through). 

At the Geodetic Institute bberehas_been developed a softwa- 
re package for gravrty ‘^‘d ^odelling^ ®AVSOFT ;Tacherning and 

and have been used *9 includes programs for gravrty mo- 
Forsberg, 1986). Ihe p ,^ pnrnT) and pFT (GEOFOUR), as well 

delling using collocation GEOCOL1. and FFT £or the 

as programs for the estimation of statistical y 
gravity field ( EMPCOV , COVFIT). 

r«D &\ 7 QnFT for aeoid dctcrnunst ion 
It is our intension to use GRAVSOFT tor geo ^ There _ 

?:re a wrwan°ted°to £ 

tuation, namely where the geoi ity da ta should be good, 

and n^ “ o^^ our^a lu^ 

^“: te su^ im f e 3itu^;o: hl is f^n^e centra 1 Mediterranean 
Sea , see Fig. 1 • 

In the following we will describe the data and the result of 
the evaluation. 


96 


ORIGINAL FACt Li 
OF POOR QUALITY 



la bio Data 


Gravity data were marie available to us by D. Arabelos, Uni-' 
vers tty of Thessaloniki, in the form of data in a 0.125° grid 
digitized from the maps published by Morelli et al.(1975). Data 
in the sea area shown in Fig. 1 was used. 

Since the use of the FFT requires data to be available in a 
regular grid all the missing values were predicted using a fast 
collocation procedure implemented in the program module GEOGRID. 
On the other hand does collocation not permit the use of all the 
4194 values, since a full system of equations with this dimen- 
sion must be solved. Therefore, when using collocation for geoid 
computation, only the 0.25° grid points were used outside the 4° 
x4° inner area, where the geoid was computed. 

Cross-over adjusted SEASAT-data (Cruz and Rapp, 1982) were 
made available to us by R.H. Rapp. A local cross-over adjust- 
ment, using the data in the 4°x4° area, made the standard devia- 
tion of the cross-over values of the six used tracks decrease 
from 0.05 cm to 0.02 cm. Raw GEOSAT data were also adjusted, and 
cross errors (mainly due to data over land, see Fig. 1) were re- 
moved. The data covered a 1/2 year period, and contained there- 
fore up to 10 repeat tracks. Originally the dataset consisted of 
3096 points, which before the removal of gross errors had a 
standard deviation of 5.53, and with 97 values removed had a 
standard deviation of 2.87 m. The result of a cross-over adjust- 
ment with only bias gave a standard deviation of the cross-over 
differences of 0.05 m compared to 5.29 m before the adjustment. 


3. Gravimetric Geoid Computations 


First the contribution of the spherical harmonic expansion 
GPM2 was subtracted. Using these "reduced" values, empirical 
auto- and cross -covariance functions were estimated by EMPCOV, 
using the gravity and the GEOSAT data, regarded as geoid 
heights. An analytic expression for the covariance function was 
then determined using COVFIT ( Knudsen , 1 987 ) , 


C< ip ) 


N 

a • I e . P . ( cosip ) 
1-0 


oo 

+ z 

i = N + l 


A ( i -1 ) 

( i-2) ( i + 4) 


R 2i + 4 

< — ) P • ( COSl/J ) . 

R l 


Here is the spherical distance between two gravity anomaly va- 
lues (at the sea surface), e^ the error degree-variances of 
GPM2 , a, A scale factors and Rg the radius of the so-called 
Bjerhammar sphere. R is the mean radius of the Earth, and Py are 
the Legendre polynomials. Values of N=12G, a=0.88, A=444 and R- 
Rg=3.75 km was found to give a nearly perfect agreement between 
the analytic expression and the empirical auto- and cross 
covariances . 


The gravity data were then used to compute geoid heights for 
the 4°x4° area. The use of collocation look more than 10 times 
as long time as the use of FFT. A comparison with the altimeter 
measurements were then made, and the results are given in Table 
1. In Fig. 2 are shown the FFT, collocation and GEOSAT heights 
along the longest track in the open sea. 
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Table 1. Comparison of If FT and collocation gravimetric yeoids 
with SEASAT and GEOSAT adjusted altimeter heights. 




Mean 

m 

Standard Dev. 
m 

GEOSAT-data 

with GPM2 subtracted 

-1.24 

0.62 

Difference 

GEOSAT- FFT geoid 

-2.18 

0.15 

Difference 

GEOSAT-Col location geoid 

-1.37 

0.15 

Difference 

SEASAT-FFT geoid 

-0.54 

0 .20 

Difference 

FFT-Col location geoid 

-0.82 

0 .20 


The difference between the FFT and collocation geoid heights 
are shown in Fig. 3. The large mean difference and standard de- 
viation may be caused by the way the two methods accounts for 
the long-wavelength information. Also the standard deviation of 
the differences is surprisingly large, considering that both 
methods agree so well with the GEOSAT data. 

A detailed analysis of the differences between the GEOSAT 
heights and the gravimetric geoid heights along the individual 
tracks, see Fig. 4, showed that altimeter data close to the 
coast <<50 km distance) have a larger variation than points at 
the open sea. This indicates a possible coastal current, the 
existence of which must be verified. 

4. Conclusion 

The result of the investigation shows (as expected) that the 
GEOSAT data in this area are slightly superior to the SEASAT da- 
ta. Also, considering the error in the altimeter data, we have 
demonstrated that it is possible to compute a regional, relative 
geoid, at the decimeter level, using the GRAVSOFT programs. It 
is obvious, that FFT should be used if the data configuration 
and quality permits it. Otherwise collocation should be used, 
since it puts few requirements on the data configuration, and 
also makes it possible to include the adjusted altimeter data as 
observations. The quite large differences between FFT and collo- 
cation must be further studied. 
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supported by European Community contract no. St2J— 39-DK { EDB ) . 
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fig. 1. Geo j d heights (contour interval 1 nl 
derived fron GPH2 with the used GEOSAT tracks. 




, collocatio. 
geoid 


2. Adjusted GEOSAT heights and corresponding 
FFT and collocation geoid heights along the longest 
track. The contribution fron GPH2 is subtracted. 


^ 2 *>■ 



Fig. 3. Difference between geoid heights 
computed by FFT and collocation. Contour 
interval 0.1 m. 



Fig. 4. Sea-surface heights computed as the difference 
between the GEOSAT heights and the gravimetric geoid 
heights computed by FFT. 
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